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Many molecular magnetic clusters, such as single-molecule magnets, are characterized by spin 
ground states with defined total spin S exhibiting zero-field-splittings. In this work, the spectro- 
scopic intensities of the transitions within the ground-state multiplet are analyzed. In particular, 
the efTects of a mixing with higher-lying spin multiplets, which is produced by anisotropic inter- 
actions and is neglected in the standard single-spin description, are investigated systematically for 
the two experimental techniques of inelastic neutron scattering (INS) and electron paramagnetic 
resonance (EPR), with emphasis on the former technique. The spectroscopic transition intensities 
are calculated analytically by constructing corresponding effective spin operators perturbationally 
up to second order and consequently using irreducible tensor operator techniques. Three main 
effects of spin mixing are observed. Firstly, a pronounced dependence of the INS intensities on 
the momentum transfer Q, with a typical oscillatory behavior, emerges in first order, signaling the 
many-spin nature of the wave functions in exchange-coupled clusters. Secondly, as compared to 
the results of a first-order calculation, the intensities of the transitions within the spin multiplet 
are affected differently by spin mixing. This allows one, thirdly, to differentiate the higher-order 
contributions to the cluster magnetic anisotropy which come from the single-ion ligand-field terms 
and spin mixing, respectively. The analytical results are illustrated by means of the three examples 
of an antiferromagnetic heteronuclear dimer, the Mn-[3 x 3] grid molecule, and the single- molecule 
magnet Mni2. 

PACS numbers: 33.15.Kr, 71.70.-d, 75.10.Jm 



I. INTRODUCTION 

Since the discovery of magnetic hysteresis of pure 
molecular origin in the molecule Mni2 about a decade 
agOfi the field of molecular nanomagnets matured 
into one of the most attractive research areas in 
magnetism^i^ In molecular nanomagnets, a defined 
number of metal ions are linked by organic ligands such as 
to form well defined, small magnetic nanoclusters. Quite 
a number of spectacular quantum effects could be ob- 
served, like quantum tunneling of the magnetization;^ 
quantum phase interference effects^ quantum rotations 
and tunneling of the Neel Yectoi^^ quantum magneto- 
oscillationSr!^ or quantum superposition of high-spin 
states lii 

On a fundamental level, the magnetism in these clus- 
ters is described by a spin Hamiltonian, which in- 
cludes the isotropic Heisenberg exchange interactions, 
the ligand-field interactions, the Zeeman term, and fur- 
ther anisotropic effects. In the majority of cases the sit- 
uation is such, that the isotropic Heisenberg exchange 
is the most dominating term, and the so-called strong- 
exchange limit, in which the non-Heisenberg terms are 
treated in first-order perturbation theory, emerged as a 
standard approach to analyze their magnetismiiSiiiiiii^ 
More precisely, a situation is considered where the dom- 
inating Heisenberg exchange results in an energy spec- 
trum whose interesting part consists of well separated 
spin multiplets - each classified by the total-spin quan- 
tum number S - which each exhibit smaller zero-field- 
splittings due to the anisotropic terms in the spin Hamil- 



tonian. A perturbational approach is then obvious. 

In the recent years, however, it became increasingly 
clear that the strong-exchange limit misses important 
aspects and effects. For instance, it has been noted 
early on that the experimentally observed tunneling split- 
ting in single-molecule magnets can be off by orders of 
magnitude from that calculated from the spin Hamilto- 
nian describing their ground-state spin multiplets, be- 
cause higher-order spin terms, most of which are exper- 
imentally inaccessible and thus disregarded, have a cru- 
cial influence liiSiii The higher-order spin terms origi- 
nate from either microscopic anisotropy terms such as 
single-ion ligand-field terms, or from a mixing between 
the different spin multiplets, which arises in second and 
higher order perturbation theory (this type of mixing has 
been denoted e.g. as S mixing in Ref. Il8|) . The origin 
of the magnetic anisotropy of the cluster is an obviously 
important question, and has been the subject of many 
studies.^^d'^i^°i^^i^^i^'^i^'^i^^ As a result of them it turned 
out, that mixing between spin multiplets can be an im- 
portant source of cluster anisotropy - and for its under- 
standing one needs to go beyond first-order perturbation 
theory. Spin mixing has been observed recently also to 
lead to another interesting effect. In a molecular Mn- 
[3 X 3] grid striking quantum magneto-oscillation were 
observed in the field dependence of the torque at low 
temperatures.^" These oscillations were shown to be di- 
rectly related to a mixing of the levels at field-induced 
level crossings, and as such may be considered as the 
most clear manifestation of mixing of spin multiplets yet 
observed. 

All these effects are intrinsically related to the fact. 
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that the energy spectrum consists of more than one spin 
multiplet, i.e., that the total spin of each spin multiplct 
arises as a result of the competing exchange interactions 
in a many-spin system. The energy levels of the indi- 
vidual spin multiplets are efficiently described by intro- 
ducing effective spin Hamiltonians for each of the spin 
multiplets of interest, which in first-order are naturally 
obtained by standard projection techniques. Here, each 
term in the microscopic spin Hamiltonian corresponds 
to a similar term in the effective spin Hamiltonian. For 
instance, a single-ion anisotropy operator such as Sf^ pro- 
duces a term a{i)S'^, where i indexes the spin centers in 
the cluster and a{i) is a projection coefficient. This ap- 
proach apparently cannot grasp the energy shifts due to 
the admixture of levels from the other spin multiplets. 
This effect leads to small changes in the coefficients a{i) 
and, more importantly, to new spin terms of higher-order. 
For example, it is intuitively clear that a microscopic spin 
term Sf^ adds a term proportional to to the effec- 
tive spin Hamiltonian, since in second-order perturbation 
theory elements like {0\SfJn){n\SfjO)/A appear, which 
behave like within the space of a specific spin multi- 
plet. This mechanism has been known for long and sev- 
eral methods have been developed to extend the effective 
spin Hamiltonian concept beyond first-order. Recently, 
Liviotti et al. presented a rather appealing strategy to 
incorporate the spin-mixing effects in the effective Hamil- 
tonian up to second-order in perturbation theoryii^ 

By construction, the effective spin Hamiltonian is de- 
signed to reproduce the energy levels of a spin multi- 
plet correctly, and thus allows to describe thermody- 
namic measurements or the transition energies in spec- 
troscopic work with high accuracy. However, this alone 
is not sufficient to ensure that also the transition matrix 
elements, which determine spectroscopically observed in- 
tensities, are reproduced correctly. To see this, let us 
consider a specific spin multiplet. Due to spin mixing, 
a small amount of the wave functions from the other 
spin multiplets is mixed in, and symbolically we write 
|0) = |0) -I- e|n), where e is a measure for the mixing (|0) 
and |n) denote unperturbed wave functions). Equiva- 
lently one may write |0) = T|0), with the transformation 
T chosen appropriately. Instead of calculating the ener- 
gies from (0|-ff |0), it is the idea of the effective Hamilto- 
nian concept to calculate it from (0|Ttffr|0) = (0|fl"|0), 

with the effective Hamiltonian H = T^HT (a short but 
excellent review of the effective Hamiltonian concept is 
found in Ref. l27j) . One may express it this way: physi- 
cally the wave functions are modified by the mixing, but 
the effect with respect to the energies can be described 
within a framework, which keeps the Hilbert space fixed 
but introduces^ a modified Hamiltonian so as to ensure 
(0|/f|0) = (0|/f|0). Since the enlargement of the space 
which is required to account for the mixing cannot be ac- 
commodated by the zero-order wave functions, it is intu- 
itively not surprising that the Hamilton operator in turn 
has to be adapted, or extended, respectively. However, 



the very same is true for any operator O if its matrix ele- 
ments are to be calculated within the restricted space of a 
spin multiplet. Therefore, similar as for the Hamiltonian, 
the operator O has to be substituted by an effective op- 
erator O, in which additional terms with no direct coun- 
terparts appear, in order to ensure {n\0\m) — {nlOlfh). 
In short, if a mixing between spin multiplets is non- 
negligible, the operators which determine the spectro- 
scopic intensities need to be adapted or extended too, in 
order to yield correct results. At this point it is worth 
mentioning that in a complete numerical diagonalization 
these issues of course do not arise, since the exact eigen 
functions obtained this way trivially include the mixing 
of different multiplets. 

In this work, we will analyze the effects of spin mixing 
on the spectroscopic intensities for transitions within a 
zero-field-split spin multiplet. This task essentially in- 
volves the construction of the corresponding effective op- 
erators, which is achieved by a procedure similar to that 
used by Liviotti et al.^^ for the construction of the effec- 
tive spin Hamiltonian (and, as a spin-off result, we also 
provide a considerable improvement of their method). 
The discussion will be conducted with the specific exam- 
ples of the single-molecule magnets, or more generally the 
spin clusters with non-compensated ground-state spin, in 
mind. For this class of compounds, the magnetic excita- 
tions within their ground-state multiplets have been in- 
vestigated intenselyii2422i22i22i^ Such studies provide de- 
tailed information on the magnetic anisotropy, which is of 
obvious importance for the understanding of their quan- 
tum properties. These systems, in fact, represent canon- 
ical examples, and the previous studies on the magnetic 
excitations within their high-spin ground states largely 
motivated this work. It will be shown, that a careful 
analysis of the transition intensities can yield crucial in- 
formation not obtainable otherwise. As a simple exam- 
ple, from the energy spectrum alone it is not possible to 
decide whether a term in the effective spin Hamilto- 
nian arises as a second-order effect from the single-ion 
terms Sf^ , or as a first-order result from single- ion terms 
Sf^. The spectroscopic intensity in principle allows to 
discriminate these two origins. 

At low temperatures, where only the ground-state mul- 
tiplet is relevant, these systems behave like a single spin. 
The corresponding effective spin Hamiltonian is known 
as single-spin model or giant-spin model. The transitions 
within the spin ground states were investigated experi- 
mentally by inelastic neutron scattering (INS) and elec- 
tron paramagnetic resonance (EPR). In both techniques, 
the measured intensities are related to the transition 
matrix elements of the local spin operators, {n\Sia\m) 
(a — X, y, z). In analogy to the spin density map, which 
corresponds to {n\Sia\n), we will call the set of matrix 
elements {n\Sia\m) the spin transition map. Clearly, the 
topology of the spin transition map is different for each 
transition. The effect of the spin mixing on the ground- 
state intensities may then be understood intuitively as 
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the admixing of spin transition maps with topologies dif- 
ferent from that of the ground state. 

The discussion will focus mainly on the situation en- 
countered in INS, as it is, when affordable, the technique 
of choice since the analysis is not affected by the com- 
plexities introduced by an applied magnetic field; the 
discussion for EPR will be limited to some important 
general issues^ INS is also of special interest since it 
allows to probe the local transition matrix elements via 
the dependence of the scattering intensity on momentum 
transfer (Q dependence), i.e., one may say that INS al- 
lows to probe the topology of the spin transition map, 
while EPR relates to an average. In spin clusters, the Q 
dependence exhibits typical oscillations due to so-called 
interference factorS)2^i2^ which are absent for mononu- 
clear compounds and are thus intrinsically a many-spin 
effect. However, the INS transitions within the ground- 
state multiplets of high-spin clusters have been analyzed 
so far as if these systems were true single spins, neglect- 
ing the Q dependence (the only exception we are aware 
of is Ref. Is^). This approach, denoted as single-spin ap- 
proach henceforth, leads to some subtleties, which will 
be addressed. 

The manuscript is organized as follows. First, in the 
next section, the single-spin approach is surveyed as this 
sets the context for the remainder of the work. Then, in 
section IIIII the problem is carefully stated and the no- 
tation set up. Also, some important relations are given, 
in particular the very useful eq. (|18|l . which is key to 
this work. The formalism developed for the calculation 
of spin-mixing effects on the transition intensities turned 
out to be very efficient also for the construction of the 
effective spin Hamiltonian, which is considered in sec- 
tion IIVI where also a general result is derived. In sec- 
tion then, which represents the core of this work, the 
effects of spin-mixing on the INS transitions are discussed 
theoretically and by means of some examples. Section lvTI 
is devoted to a discussion of some effects of mixing in 
EPR. Section IVIII finally presents a conclusion. Some 
supporting information, helpful for the understanding of 
the formalism and for actual calculations, is provided 
in the appendix. Readers who are mainly interested in 
a qualitative understanding of the effects may skip the 
technical subsections marked by an asterisk. 



II. SURVEY OF THE SINGLE-SPIN 
APPROACH FOR ANALYZING INS DATA OF 
HIGH-SPIN CLUSTERS 

The INS transitions within the ground-state spin mul- 
tiplet of high-spin clusters were analyzed so far using a 
formula derived for mononuclear compounds (with the 
exception of Ref. f36l) . For a mononuclear paramagnet, 
the transition intensity from state |n) to |m), Inm, is pro- 
portional to the matrix element |(n|S'^|m)p, where is 
the spin component perpendicular to the scattering vec- 
tor Q = k — k'j^ For a powder sample, after averaging 



over all orientations of Q, one obtains 

/„,„ cxp„ (2|(n|^,|m)|2 + \{n\S+\m)\^ + |(n|5_|m)|2) 

(1) 

(p„ is the Boltzmann factor for the state n)?^-'^'^ 

The use of this formula was certainly motivated by the 
insight, that high-spin clusters in many respects behave 
at low temperatures like large, single spins. However, us- 
ing eq. , which is the single-spin approach, can produce 
substantially erroneous results for the calculated INS in- 
tensities. 

It is a celebrated feature of the INS intensity to pro- 
vide information about the wave functions (for a recent 
discussion see Ref. |3^ and the references cited therein) . 
It is thus intuitively clear, that the INS intensity "feels" 
the differences between a single-spin and a many-spin 
wave function. In fact, for a polynuclear spin cluster 
the INS intensity depends on the local transition ma- 
trix elements and the geometrical arrangement of the 
spin centers in the cluster via a factor exp[iQ(Ri — 
T{.j)]{n\Sia\m){'m\Sj[3\n) (Ri denotes the position vector 
of the ith spin center). This gives rise to characteristic in- 
terference factors in the INS cross section of powder sam- 
ples, and hence characteristic Q dependencies of the INS 
intensity, which are typically oscillatoryi^^iS^ For a single 
spin, in contrast, such interference effects are absent, and 
the INS intensity is independent of Q, see eq. (QJ. Thus, 
as mentioned already in the introduction, the interfer- 
ence effects visible in the Q dependencies are a direct 
fingerprint of the many-spin nature of the wave func- 
tions in polynuclear clusters. These effects were indeed 
observed for, e.g., the single-molecule magnet Mni2r& 
showing that the deficit of the single-spin approach to 
reproduce the Q dependence is experimentally relevant. 

In a first approximation, the Q dependence is the same 
for all transitions within a spin multiplet. In fact, in sec- 
tion we will show that this is strictly correct up to 
first-order in perturbation theory. Thus, if one is analyz- 
ing the transition intensities at constant Qi the relative 
intensities of the transitions as calculated from eq. (Q) 
are reasonably accurate. However, due to obvious ad- 
vantages, the spectra are often recorded on time-of-flight 
(TOF) spectrometers, and the integrated spectra, which 
are obtained by summing the data of all detector banks, 
are used in the analysis. In that case a rather trivial but 
important problem occurs, which we mention here since 
it seems not to have been considered in the previous anal- 
yses: 

In the integrated TOF spectra, the intensity of each 
peak corresponds to 7*'^ = ^^^"^^ g{Q)In^{Q)dQM 
Thus, even if the INS intensities are independent on Q, 
the TOF intensities are not simply proportional to the 
squared matrix elements and the Boltzmann factor, but 
also scale with G{E) = ^^^Z diQW- ^i^ce G{E), how- 
ever, depends only on instrumental parameters, it is 
easy enough to account for that (though this has been 
disregarded in previous works using the single-spin ap- 
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proach). But this procedure is not working if the Q 
dependence exhibits pronounced maxima and minima, 
as it is often the case. Then, for one peak Qmm and 
Qmax may be such as to bracket just the top of a maxi- 
mum, while for another peak they may bracket also the 
nearby minima - with corresponding effects on the inte- 
grated peak intensities. In conclusion, when using the 
single-spin approach for the analysis of integrated TOF 
spectra one should be aware of a possibility of signifi- 
cant deviations between experimental and calculated in- 
tensities. As a side-effect, the TOF intensities of the 
Stokes and anti-Stokes peaks will not obey the relation- 
ship /(— w) = exjp{fiw/kBT)I[Lij). The situation can be 
improved considerably if a rough approximation for the 
Q dependence is available, e.g., as extracted from the 
TOF data itself, or if approximations for the wave func- 
tions are conceivable. In the remainder of this work, we 
will disregard these issues, but focus entirely on the ef- 
fects of the many-spin nature of the wave functions on 
the intensities. 

Finally it is mentioned that the above considerations 
are rather specific to INS. For instance, the EPR signal 
does not provide information about the spatial proper- 
ties of the involved wave functions as it uses spatially 
homogeneous magnetic fields Accordingly, the single- 
spin model works much better for the EPR intensity. In 
the language of perturbation theory, the Q dependence of 
the INS intensity is a first-order effect of the many-spin 
nature of the wave functions, while the EPR intensity is 
affected only in second-order {vide infra). 

III. GENERAL RELATIONS 

The microscopic spin Hamiltonian describing a mag- 
netic spin cluster of N spin centers, where the ith center 
has spin length Si, is 

H ^ — ^ Jij Si • Sj + ^ Si • n,;j • Sj 

-f ^ S. • D, . S, + ^ ^ B^(i)6l\i) 

% i m 

+/.B5ZS^-g. -B. (2) 

i 

The first term represents the isotropic exchange inter- 
actions, the second term describes anisotropic and anti- 
symmetric exchange and dipole-dipole interactions, the 
third and fourth terms are single-ion ligand field terms, 
and the last one is the Zeeman interaction. This expres- 
sion is quite general; from the usually considered terms 
only the higher-order exchange terms like biquadratic ex- 
change were excluded (and this only for convenience). 
The tensors are traceless and symmetric, 0^{i) are 
Stevens operators acting on the ith spin,"41' and is the 
Bohr magneton. 

In this work we are considering systems in which the 
isotropic exchange interaction dominates the other terms. 



Then the basis states of the Hilbert space are conve- 
niently chosen as the spin states \kSM), where k ab- 
breviates the intermediate spin values according to the 
chosen coupling scheme. S and M denote the total spin 
and magnetic quantum number, respectively. Working 
in this basis, the irreducible tensor operator (ITO) tech- 
nique provides a powerful tool for the calculation of ma- 
trix elementsiiSi^ In terms of ITOs, the above Hamilto- 
nian writes 

+ Erf *(D,)ff w + ErW*(S4W)TW(0 

i,q i,q 

+ J^T('>{h,)f!,'H^). (3) 

Here, hi = hbSi ' B was introduced. Tq''\Y) denotes 
the qth component of the irreducible spherical tensor 
of rank k related to the Cartesian tensor V [it is pro- 
portional to the spherical harmonics Ykq^V), and obeys 
T^''^* = (-1)it[^^ ]. f^''\i) denotes an irreducible tensor 
operator which acts in the spin space of the spin center 
i, and Tq^\ij) denotes an ITO which is composed of two 
single-spin ITOs acting in the spin spaces of the spins i 
and J with i ^ j, respectively^^ In equation each 
term of eq. is replaced by a corresponding spherical 
scalar product, expressing the fact that a Hamiltonian 
transforms under rotations in space as a scalar (see ap- 
pendix for more details concerning the Cartesian and 
spherical representations of a spin Hamiltonian). 

The second and third terms in eq. (0) or eq. (O, re- 
spectively, can be combined into J2ij S« ' ^ij ' ~ 
^^^.^rf ^*(A,,)ff )(zj) if one defines A,j = n,y(l-<5,y) + 
DiSij, and adopts the convention that Tq^\ii) means 
Tq^\i). This will simplify calculations considerably as 
the bilinear anisotropy terms in H can be treated simul- 
taneously this way. 

Many results in this work depend only on the trans- 
formation properties of the ITOs, and not on ij or i, i.e., 
on which spin centers are involved. We often simplify 
notation by writing Tq^\s), where s stands for either ij 
or i (and similarly in sums like X^s)- meaning is 

obvious from the context, we also write just Tq ' . 

The perturbative approach which is adopted to calcu- 
late energies and transition matrix elements is as follows. 
The unperturbed Hamiltonian consists of the isotropic 
exchange interactions. Ha = — J^ij Jij^i ■ Sj. The cor- 
responding zero-order eigenstates are written as \tSM), 
where r denotes additional quantum numbers. 

The aim is now the following: Within one of the (non- 
degenerate) spin multiplets specified by tS, the micro- 
scopic spin Hamiltonian H shall be replaced by an ef- 
fective spin Hamiltonian operating in a single-spin space 
\SM). This is achieved in two steps. First one deter- 
mines the effective Hamiltonian H acting in the subspace 
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\tSM) by a perturbation technique, and then one ex- 
presses H by single-spin operators using ITO techniques. 
Similarly, the spin operators Sia shall be replaced by ef- 
fective spin operators Sia expressed by single-spin oper- 
ators Sa- In the following subsections we provide some 
background as needed in this work, and present the very 
useful eq. H18|l (readers not interested in details may skip 
these technical sections). 



A. Effective Hamiltonian, effective operators* 

In the spectral representation, the Hamilton operator 
H becomes 



H 



\n){n\H\n){n\, 



(4) 



where \n) are the eigenfunctions of H. The effective 

Hamiltonian H describing a subset of levels of H is then 
constructed by introducing a model space §, the states 
of which will be denoted as with a one-to-one corre- 
spondence to the subset of levels of H, i.e., n<-^ n. Then 



J2\n){n\H\n){n\ 



(5) 



where Ps — introduced. Frequently, as is 

the case in this work, the model space § is a subspace 
of the Hilbert space of H, then the procedure consists of 
finding an unitary transformation with \n) — T\n): 



{n\H\n) 



\t''ht[ 



(n\H\n) 



(6) 



Several methods are available to find the transforma- 
tion T perturbativelyiSi^Si£Li2ii2ia As usual, one starts 
from H = Hq + Hi with Hi being small. In general, the 
spectrum of Hq exhibits degeneracies. The eigenstates of 
Hq are thus written as \lfi) such that Hq\1^) = Ei\lfi). 
The different methods lead to equivalent results in the 
order considered here (the results below seem to have 
been first reported by Pryce'*^). Following Schrieffer and 
Wolff,^^ or Schlichter,"''" one chooses T = e^^ , which im- 
plies H ^ H + i[H,U]-l/2[[H, U],U], and separates the 
perturbation as Hi = Ha + Hi, such that {l^jL\Ha\l' fJ-') — 
for I ^ /' and {l^i\Hb\l' f^') = for / = I'. Then, the 
matrix elements of Ht can be minimized by choosing 
Hb + i[Ho + Ha, U] ~ 0, which results in 



H^Ho + Ha + i/2[Hb,U] 



(7) 



up to second order in Hi, and/or U, respectively. The 
transformation U is obtained as 



_ 1 PiHbPi' 

i ^ Ei> - Ei 

i^i' 



(8) 



and the effective Hamiltonian acting in the subspace Pi 
as 



(9) 



with Pi = Note that PiHPr = PiHbPv for 

I ^ v. 

Similarly as for the Hamiltonian, an arbitrary oper- 
ator O has to be replaced too by an effective opera- 
tor O acting in the model space § in order to ensure 
{n\0\m) = {n\0\m). Accordingly, 



6 = PsOPs = \n){n\d\m){f 



(10) 



or, if the model space § is again a subspace of the Hilbert 
space of H, 



O = e"*^Oe 



(11) 



In contrast to the result for the effective Hamiltonian, 
the effective operator is obtained as O = O + i[0,U]. 
The difference occurs because Hb is of lower order than 
O. Inserting the above perturbation result for the trans- 
formation J7, one finally obtains 



O 




HPi,0 + OPvH 

Ev — El 



Pi- 



(12) 



In the context of this work, having written the zero- 
order eigenstates as \tSM), the index I corresponds to 
tS*, and /i to M . The operators obtained from eqs. 
or (|12|l . which are acting in the subspace \tSM) with r 
and S fixed, can be replaced, as will be discussed in the 
next subsections, by operators acting in the single-spin 
space \SM). 

In order to avoid confusion which spin space is meant, 
we apply the following conventions: Whenever both 
quantum numbers r and S appear in a matrix or reduced 
matrix element, it should be understood as a many-spin 
matrix element, while if only S appears, a single-spin 
matrix element is meant. If not already unequivocally 
clarified by this, Tq'^\s) indicates a many-spin ITO act- 
ing in the space \tSM), and Tq''\s) a single-spin ITO 
acting in the space \SM). 



B. Expansion of a spin operator into ITOs* 

The aim is to find the expansion of an operator 0(5*), 
which is a function of spin operators, into irreducible 
tensor operators: 



kq 



(13) 
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The coefBcients can be calculated by considering the ma- 
trix elements {SM\0\SM') within a spin multiplet. With 
the help of the Wigner-Eckardt theorem and the orthog- 
onality relation 



E 



Jl J2 J3 

mi m2 TO3 



Ji J2 J3 

mi 7712 



where (...) denotes a Wigner-3j symbol, one directly ob- 
tains the following useful result: 



Ckq 



Em; 



M-S 



MM' 



S 
-M 



k S 
q M' 



2fc+ 1 
~ (5||f(fe)||5) 

{SM\d\SM'). (14) 



Here, (S'H'f^'^^llS') denotes the (single-spin) reduced 
matrix element. If the operator O is acting in 
a many-spin space, 0{S) should be interpreted as 
E MM' I > {rSM \d\TSM'){SM'\. 



C. Spin operator equivalents* 

The procedure to find the spin operator equivalents of 
the first-order terms PtsOPtS is well known, -'^^ and is 
reproduced here for convenience. One first expresses O 
by (many-spin) ITOs, similar as in eq. (|13() . Then the 
Wigner-Eckardt theorem is exploited resulting in 

Prsf^''\s)PrS = rfe(s)fW(5), 



rS\\f^^\s)\\TS) 

{S\\f(^)\\S) 



(15) 



with the projection coefficients Tk{s) (in case of ambigu- 
ities, Ffc will be called a f^'-order projection coefficient). 

To find the spin operator equivalents for expressions of 
the form PtsOiPt'S'02PtS, which appear in second or- 
der, one again first decomposes the operators Oi and O2 
into ITOs. This leaves one with the evaluation of terms 
alike PrsT^j^^^ Pr'S'T^i^^^ PrS- Though they could be re- 
placed directly by effective single-spin operators (which, 
however, are not easily found), it is more convenient to 

expand PrsTq^^'' Pt' S'Tq^^'' PtS in turn into single-spin 
ITOs because, as will be shown below, the respective ex- 
pansion coefhcients can be calculated very easily. 

The matrix elements oi Prsf^t"-^ Pr'S'T^'^''^ PrS within 
the spin multiplet labelled by tS are 

^{TSM\f^'^^^\T'S'M') {T'S'M'\f^^^^\TSM") = 

M' 

(T5||f(^i)||r'5')(T'5'||f('=^)||T5)^(-l)^-^'^+^'-*^'x 

M' 



S ki S' 
-M qi M' 



S' k2 S 
-M' q2 M" 



Insertion of this expression into eq. (|14() in order to get 
the expansion coefficients requires the evaluation of 



-M' I S k S 
-M q M" 



E (-1)^' 

MM'M" 

S ki S' \ f S' k2 S 
-M qi M' I \ -M' q2 M" ' ' 



(16) 



which, by using some symmetry properties of the Wigner- 
3j symbols and (-1)"^^' = {-ly^^'+^t+^t" -11+12 ^ can be 
transformed into 



(-1) 



S'-gi+g2 



E(_T\M'+M+M" I S' S ki 
^ > \ M' ~M qi 



MM'M" 

s s 



M -M" -q \ M" -M' q2 



With the help of a relation between the Wigner-3j and 
6j symbols. 



E (-1 

M1M2M3 



,.7i + ,/2 + ./3+Mi+M2+Af3 ( J'i 33 

Ml -M2 ms 



J2 J3 ji \ ( J3 h 
M2 -M3 mi J [ M3 -Ml TO2 



Jl J2 J3 
mi 7712 fTT-S 



where {...} denotes a Wigner-6j symbol, eq. ((TB|l assumes 
the simple expression 





h 


h \ 




J2 


J3/ 



f k k2 


ki \ 




k2 


ki \ 


V -1 12 


qi J 


is' 


S 


S j 



Using (tSIIT'^Wt'S') = (-l)^'-^^(r'5'||f('=)t||TS')* and 
introducing the coefficients 

Ck.k.k = (-l)^+^'(2fc-t-l)(5||fW||5)| ^, 



rfcifc(s) 



ki k2 k 

-qi -92 q 

(r5||f(^-i)(s)||r'5'> 



(5|lfW||5) 



(17) 



where we call Tj^-^kis) a 2"''-order projection coefficient, 
a very useful formula is finally obtained: 



Prsf^'' Hsi)Pr'S' f^^'^ {s2)PrS^ 
Y.^kA-^l)nA'2)4lq^qT(''HS). 



(18) 



Here, the sum runs only over k because of the restric- 
tion q = qi + q2- The similar structure of eq. (|18|l and 
eq. ifT^ is satisfying. 
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The coefScients d'^lq^q do not depend on the actual 
wave functions involved and thus can be calculated, and 
tabulated, once and for all. The potentially most needed 
coefficients are listed in TableQ] For calculations it is use- 
ful to realize that with respect to the first two " columns" 
they have the same symmetry properties as the Wigner- 
3j symbols, e.g., d^^^^^ = {-l)'''+''^-^''d'^l^^^ . It is also 
useful to note that for a given ki the 2"''-order projec- 
tion coefficients for different k are trivially related to each 
other as they all refer to the same reduced matrix ele- 
ment, e.g., ri,fc(s) = (5||f (i)||^)/(5||f W||5)ri,i(s). 



IV. EFFECTIVE SPIN HAMILTONIAN* 

Using the results of the previous section, a general 
equation can be derived for the second-order effective 
spin Hamiltonian. A general procedure, which in princi- 
pal leads to the same results as ours, has been described 
recently by Liviotti et al.^^ It consists of introducing new 
spin operators -^^fe(e/c)' '^hich in turn can be written 
as a combination of Stevens operators. However, no sim- 
ple and straightforward procedure is available to obtain 
the M operators. Furthermore, as the analysis of exper- 
imental data is generally done in terms of the Stevens 
operators (or equivalents), one almost always wants to 
have the ^'^k{e/c) operators expressed by Stevens op- 
erators anyhow. We believe that our approach has some 
advantages: it produces the effective spin Hamiltonian 
directly in terms of ITOs (to which the Stevens oper- 
ators are simply proportional), and, more importantly, 
eq. H18|l of the previous section provides a pretty easy 
way to find this decomposition. 

The calculation of the first-order contribution to the 
effective spin Hamiltonian via eq. (|15|l is standard and 
needs no further comment,-- The calculation of the 
second-order contributions is enabled by noting that the 
perturbative part Hi of the microscopic spin Hamilto- 
nian, eq. consists of a sum of spherical scalar prod- 
ucts. Hi may thus be decomposed as 



Hi — Hr, 

r 

= EE^.^^*(V,jT(^)(s.), (19) 

where r (and u) indexes the individual terms in 
Hi. The second-order contribution produced by one 
higher-lying spin multiplet is thus proportional to 
J2r,u PrsHrPr' s' HuPtS, See eq. ©. Inserting the indi- 
vidual terms of this sum in eq. (|18|l leads to expressions 
like 



Erfc.fc(Sr-)rfc^fc(s„) Cfc^fc„fc E ( 
k q^qu \ 



-Qr -Qu q J 

X r!t^ ( v,,_ )t['-^ ( )rf ) (s) . (20) 



The sum over qr and is nothing else than the tensor 
product of two irreducible spherical tensors. With that, 
the effective spin Hamiltonian up to second order may be 
conveniently written as 



kq 



with 

"-qr 



"■q 



(21) 



(-1) 



kr-k 



^ (V, J T^*^") (V, J] W . 



1 



(22) 

Here, only one higher-lying spin multiplet with a gap A 
has been considered, but the generalization to more spin 
multiplets is obvious [only multiplets with jS* — S"| = 

0, . . . , min{kr, k^) contribute to h^^^^]. The similarity of 
the two terms in eq. ^21\ with spherical scalar products is 
satisfying to note. Unfortunately, such an elegant result 
cannot be obtained for the effective spin operators Sia- 
Here terms like PtsHtPt' S'Tq^\s)PrS appear in which 
only one sum over a component, q^, is involved. 

These equations shall be applied now to the case where 
the small part of the microscopic Hamiltonian consists of 
only the bilinear anisotropic term, i.e.. Hi = S^- A^j • 

Sj , were the tensors A^ j are assumed to be diagonal and 
traceless. This situation is the one most often encoun- 
tered in practice [in this case eqs. (|21|l and (|22|1 simplify: 
kr and k^ are replaced by 2, Sr and s„ by s, and remain- 
ing indices r,u can be dropped]. The tensors A<, can be 
parameterized as usual by Dg and Eg parameters, and 
the non-zero components of the corresponding spherical 

n(2) 



andT2^"^(A,) 



tensors become Tq (A^) 

(see also appendix^. After a simple, but lengthy cal- 
culation the effective Hamiltonian eq. I|21|) is obtained 
as 

h = D[Sl-l/?,S{S + E{Sl- sD 

O/SO , d2A2 , d4A4 



D = D 



-BlOl + BiOi + BlOl, 
- ^ d^^[{D,f + i{E,n, 



E = E 



V3A 
V2 



^^df^iD^E* + D*^E^), 



1 



6V70A 



dlll[2{D,f + {E,f% 



^ dllliD^El+DlE^l 



Bi = - 



3V14A 



6V2A 



^^000 (-^4)^- 



(23) 



Here, we have expressed ITOs by the conventional op- 
erators using Table to put the effective Hamiltonian 
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TABLE I: The coefficients dooo^'' discussed in ttie text for the potentially most need cases. Independent on S' — S is dgoo = 1/5. 



5" 






"000 


^00A 

"000 


j121 
"000 


"000 


5"- 


2 




+ l)(2S + 3) 


|^(5+l)(2S + 3)(5 + 2)(2S + 5) 


^V(5-i)(s + i) 




S- 


1 


^is- 


-5)(2S + 3) 


-4#(25 - 3)(25 + 3)(5 + 2)(2S + 5) 


-4oo^(2S + 3)(5 + 2) 


S 






-3)(25' + 5) 


^{23 - 3)(5 - 1)(5 + 2)(25 + 5) 


iv'(25-l)(25 + 3) 


dooo^(5-l)(5 + 2) 


S + 


1 


^{2S 


-l)(S + 6) 


-4#(25 - 3)(25' - 1)(5 - 1)(25' + 5) 


^VSiS + 2) 


-^000^(5-1) (25-1) 


S + 2 


21 


(25-1)5' 


|^(2S-3)(25-l)(5-l)S 







in the standard form (the constant fc = term has been 
dropped). Also, we have introduced the definition 



\-j2 



(gi?y)^Tf) *(R, 



(24) 



where the index k might be dropped, which means then 
D = DsT2{s) (analogous definitions were used for 
the Ek terms). It should be recalled that s stands for the 
indices ij or i, and that the sum over s in the present 
case actually means J2s ~ Again, only one 

higher-lying spin multiplet with a gap A has been con- 
sidered. For the special case of S" = — 1, cq. was 
derived already in Ref . . 

As expected, the mixing of spin multiplets due to the 
anisotropy leads to a renormalization of the first-order 
zero-field-splitting parameters D and E, and to the ap- 
pearance of fourth-order spin operators with strengths 
B^, Bl, and Bf. In case of metal centers with spins 
larger than 3/2, the single-ion fourth-order spin terms 
in the microscopic Hamiltonian H also contribute to the 
fourth-order terms in the effective spin Hamiltonian in 
first order. This adds terms B^' = T4{i)B^'{i) to the 
Bf's in eq. (ESJ. 



V. INELASTIC NEUTRON SCATTERING 

In this section, we analyze the effects of spin mixing 
on the INS intensity of a powder sample in zero mag- 
netic field. A discussion of the situation encountered in 
EPR, which additionally requires the consideration of a 
magnetic field, is deferred to the next section. The INS 
intensity in a spin cluster is given by the following cross 
section 



dQduj 



(25) 

were C(Q,T) = {-ie^ /m^c^){k' /k)eiq^[-2W[Q,T)], (i = 
l/iksT), and Z(T) is the partition functiomSJ-S The 
sum runs over all transitions \n) |m). For a powder 
sample in zero magnetic field Inm is given by 

InrniQ) = ^ J^* (Q)^; (Q){ | JO ) • 



t(i)(s,)®tW(Sj) } 



(2) 



(26) 



In this equation, Fi{Q) is the magnetic form factor of 
the ith spin center, jk is the spherical Bessel function of 
order k, and Ry = Rj is the distance vector between 
the ith and jth ion [explicit expressions of eq. (|26|l are 
provided in appendix [C]. The ordered products SiaSjp, 
which appear in eq. H26|l . stand for {n\Sio]rn){rn\Sji3\n), 
where \n) denotes an eigenstate of the Hamiltonian H . 

For the evaluation of the INS intensity of the transi- 
tions within a spin multiplet, one inserts {n\Sia\nn) = 

{n\Sia\m) into eq. H26|) . determines the effective spin op- 
erators Sia , and then substitutes them by single-spin op- 
erators Sa- Unfortunately, no compact result is avail- 
able for the effective spin operators, and accordingly also 
not for the INS intensity. The calculations are pretty 
straightforward, but lengthy, as are the results. Thus, in 
the next subsection (which may be skipped by readers 
not interested in the details), we limit our calculation of 

the Sia to the case of a bilinear anisotropy described by 
diagonal and traceless tensors Aij. It is, however, general 
enough to account for most situations met in practice. In 
order to gain a feeling for the effects to be expected, in 
particular as compared to the single-spin approach, we 
then - in a second subsection - consider shortly the first- 
order effects and then specialize to the case of a uniaxial 
anisotropy, for which we work out the INS intensity up 
to second order analytically. In a third subsection we 
finally discuss three examples exhibiting uniaxial mag- 
netic anisotropy. Since in this work we are interested in 
transitions within an individual spin multiplet, an uniax- 
ial system is the minimal model of relevance (concerning 
transitions between different spin multiplets, we mention 
e.g. Ref.m. 



A. Effective spin operators* 

The effective spin operators Sia within a tS subspace 
shall be calculated for a system described hy H — Hq + 
Hi with Hi = '^ij Si ■ Aij ■ Sj , where the Aij are assumed 
to be diagonal and traceless as in section Hvl According 
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to eq. (|12() . the calculation of the second-order contri- 
bution involves the terms PrsTq^\s)PT'S'Tql\i)PTS + 
Prsf^l\i)Pr'S'f^?{s)PrS- In employing eq. the 
values of k are restricted to k = 1,2,3. Because in our 
case of an isotropic zero-order Hamiltonian Hq the re- 
duced matrix elements may be chosen as real, and be- 



cause of d'^lg^q = i~^)''~^^d.l2qiq' the fc = 2 contribution 
cancels out (we nevertheless keep complex writing since 
sometimes it has advantages to work with complex wave 
functions, e.g., when exploiting spatial symmetries'^--^—). 
Re-expressing the coefficients d^'^q^q by dgoo j the effective 
spin operators finally become 



VGA 



2V15A 



[55^2 - S{S + 1) - 55^ + 2] 



HI ^= [Jx^y "r JyOx)^!i- 



d}^oDS) 




15A 



[bSl 




^dlllMi) 



15A 



2VT5A 



{Si 



(27) 



Here again, only one higher-lying spin multiplet at energy 
A was assumed (only states with jS" — S"| — 0,1 may 
contribute). Also, we have introduced the abbreviation 



(28) 



with Dk defined as in section llVl eq. H24(l (and analogous 
definitions for the Ek terms). 

The various terms appearing in Sia shall be analyzed. 
In each of the above expressions, the first terms, which 
are the first-order contributions, are constants within a 
tS subspace (or \SM) space, respectively) and equal 
for all components a = x,y, z. The second terms are 
also constants within the \SM) space, but introduce an 
anisotropy in the intensities as their signs and magni- 
tude are different for the x, y, and z components. The 
remaining terms, finally, depend via the S'a-dependent 
prefactors explicitly on the particular states involved in 
a transition, and thus lead to differences in the relative 
intensities of the transitions within a spin multiplet. 

In principle, if ions with spins larger than 3/2 are in- 
volved, there would be also contributions from the fourth- 
order single-ion terms in the microscopic spin Hamilto- 
nian. However, in contrast to the situation for the ef- 
fective spin Hamiltonian where these terms had to be 
included (though in first order was sufficient), they can 
be safely neglected here, because they appear in second 
order as the bilinear terms, but are significantly weaker 
than those. 



B. INS intensity 

We are now prepared to discuss the intensity of the INS 
transitions within a spin multiplet. First, we consider the 
scattering intensity up to first-order only, as this allows 
a general result. In this approximation, the effective spin 



operators are simply given by Sia — ^i(i)Sa [the projec- 
tion coefficient ri(i) is defined in eq. H15|) ]. which after 
insertion into the INS formula eq. H2t)|l yields 



-{-V3)CiQ)S-S 



+ E/f (Q) pi)(S)®TW(S)J , (29a) 

ftHQ) - Ei^*(Q)^;-(g)j;o(Qi?.,) 

(29b) 



xrW*(R,,)ri(z)ri(j), 



where the factor — \/3 in the first term on the r.h.s. of 



f(0) 



eq. (I29a|l compensates for Tg in the definition of /q 

[since Tl^°\Rij) = -l/VS]. The functions f^''\Q) may 
be considered as the (first-order) interference factors ap- 
propriate for the problem, as they produce a Q depen- 
dence. Appearing in first-order, they are actually the 
main source for the Q dependence of the INS intensity. 

The first-order result eq. is interesting because it 
clearly demonstrates the deficits of the single-spin ap- 
proach. Only the first term on the r.h.s. of eq. H29a|) 
is retained in the single-spin approach, but with the in- 
terference factor /q°^((5) omitted (it may be noted that 

2S • S = 2\{n\Sz\fh)\^ + \{n\S+\m)\^ + \{n\S-\m)\^). The 
effects of neglecting the Q dependence were discussed al- 
ready in section^ It should be noted that for transitions 
within a spin multiplet, the k = 2 terms in eq. H29|l are 
always of importance so that the Q dependence missed 
in the single-spin approach is not only that of /q^-* (Q) [an 
example is given below through eq. (|30|l ] . 

Equation H29|l becomes quite useful when the Fi (i) are 
known, e.g., from approximative wave functions or nu- 
merical calculations, as it then allows to calculate the INS 
intensity correctly up to first-order, which often is a great 
improvement as compared to the single-spin approach. 
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A further interesting possibility emerges. In general, the 
projection coefficients Ti{i) are affected by the exchange 
couphng topology as employed in Hq = —Jij ^ - ■ Si ■ Sj. 
An analysis of the Q dependence in principle thus pro- 
vides information about the exchange coupling constants 
in the system under consideration. 

We proceed with the calculation of the INS inten- 
sity up to second order for a system with uniaxial mag- 
netic anisotropy, i.e., a system described by the micro- 
scopic Hamiltonian H = Hq + Si ■ Aij ■ Sj, where 
the Aij are diagonal and fulfill A^x = ^yy = —^zz/'^- 
The second-order effective spin Hamiltonian is then H = 
D[Sl - 1/3S{S+ 1)] + B^O^ with parameters as given in 

eq. 12311 . Furthermore, the eigenstates of H are simply 
the \SM) states, and the calculation of matrix elements is 
elementary (in general \n) will be a linear combination of 
\SM) spin states with different magnetic quantum num- 
bers M). 

As only transitions with M M+1 (and M) 
are relevant, the a = z terms drop out of the formula for 
the INS intensity, eq. H2t)|) . One is left with considering 

Six and Siy. For a uniaxial system their structure is of the 
form Si,x/y = [ai + bi — CiF{Sz)]Sx/y, with site-dependent 
constants a.^, bi, a, and F{Sz) = 5S'^-S{S+l)-5Sz + 2, 
see eq. ^7^. The matrix element SixS-jx = {M\Bix\M + 
1)(M -I- l\Sjx\M) equals, up to first order in the gap A, 
[aiUj + aibj + biUj — aiCjF{M + 1) — CiajF(M)]SxSx, and 
similarly for SiySjy. Because of the sum J2ij in the INS 
intensity formula, and because of ^j^- aibj = J2ij biO,j 
and similarly for aiCj , one may simplify to {aiQj + 2aibj — 
a,Cj[F{M + 1) + F{M)]}SxSx. With these and similar 
results for S,ySjy, with F{M + 1) + F{M) = 2[?>NP - 
5(5+1) + 2], ~Sx~Sx=~SySy = [5(5-f l)-M(M + l)]/4, 
and the explicit expressions for the constants Oj, bi, Ci 
[eq. H27|) ]. one finally obtains 



ties of the projection coefficients Ti^kii) the relation 
53/3 = i{S\\fW\\S)/{S\\f(-^^\S))^Dj, holds, i.e., that 
up to a scaling factor fi{Q) and fsiQ) are the same. 

The function fmiQ) may be again considered as the 
appropriate interference factor, and eq. H30|l demon- 
strates all the effects which we wish to point out in this 
work. The first term in fuiQ), f{Q), corresponds ex- 
actly to the first-order contribution calculated in eq. H29() 
[/(g) = -V3/^°^(Q)-3/2/o^^^(Q)]. As expected, it leads 
to a Q dependence which is equivalent for all transitions 
within a spin multiplet, i.e., is independent on M. The 
spin mixing manifests itself in two ways. First, it adds to 
the overall Q dependence via the second term in fM{Q)- 
This contribution is independent of M like the first-order 
contribution, but smaller by a factor of 1/A. Thus, from 
an experimental point of view, it is a minor correction 
which we expect is difficult to be separated experimen- 
tally from f{Q). The second effect of spin mixing, how- 
ever, brought in by the third term, leads to different Q de- 
pendencies of the individual transitions because of its ex- 
plicit dependence on M. This effect should be amenable 
to experimental detection as a more sensitive differential 
data analysis can be applied. For instance, one could an- 
alyze the differences Imm+i/[S{S + 1) - M{AI + 1)] - 
I]i?[M'+i/[S{S + l)-M'iM' + l)]. 

The second-order projection coefficients Ti^k{j) 
are proportional to the reduced matrix elements 
(rS'llf (i)(j)||r'5') [eq. and thereby also to what we 

called the spin transition map. Thus fi{Q) and fsiQ), 
respectively, arise from the correlation of the geometrical 
structure with the spin transition map. Equation H3()|l 
confirms the intuitive expectation that a mixing with 
higher lying spin multiplets, which mixes in spin tran- 
sition maps with different topologies than that of the 
ground-state multiplet, leads to different Q dependencies 
for each transition within a ground-state multiplet. 



-Jm{Q) [5(5 + 1)-M(M + 1)], 



/m(0) = !{Q) 



MQ) 



fk{Q) = 5]i^*(g)F,(Q)[jo(QA 



V6A 

[5M^-S{S+l) + 2]h{Q), 



.h{QR^,)cU'R^^,)] ri(i)ri^fc(i). (30) 



Here, Im^^i^m — ^m^m+i ^'q (Rij) — [3{Rij^z/ Rij)'^ 
l]/2. In the definition of the terms fk{Q), the in- 
dex k might be dropped, which means to use ri(j) 
instead of Ti^kU)- The coefficients and dJ§Q are 
given m Tab. H Di and D3 are defined in eq. (|24|l . 
and the projection coefficients Ti,k{i) in eq. ((T7|l . We 
have assumed Dk and Ti^ki'i) to be real for simplic- 
ity. It is useful to note that because of the proper- 



C. Examples 

In this section, three examples of molecular spin clus- 
ters, a fictitious but instructive hetero-nuclear dimer, the 
Mn-[3 X 3] grid, and the Mni2 cluster, will be discussed 
in order to underpin the above considerations. 

As a first example a ferromagnetic spin-5/2 - spin- 1/2 
dimer with uniaxial anisotropy is considered. This ex- 
ample may seem somewhat artificial, but many hetero- 
nuclear dimers exist, and dimer clusters may evolve into 
model systems for the study of the many-spin effects on 
the INS intensity as the analysis becomes particularly 
straightforward. The microscopic spin Hamiltonian is as- 
sumed 'AS H ^ - JSi • S2 + DiSf ., with Si = 5/2 and 
^2 — 1/2. The energy spectrum is simple, it consists of 
a S* = 3 ground-state multiplet and a higher lying 5 = 2 
multiplet at energy A = 3J, which both are zero-field 
split due to the anisotropy. In the following, the values 
J = 5 K and Z?i = 1 K will be assumed. For these pa- 
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FIG. 1: Results for the ferromagnetic spin-l/2-spin-3/2 
dimer with uniaxial anisotropy considered in the text, (a) 
and (b) show the Q dependence of the interference factors 
/m for the three transitions ^ 1, 1 — > 2, and 2^3 within 
the 5* = 3 multiplet as calculated by second-order perturba- 
tion theory, eq. I30II . and numerically from the microscopic 
spin Hamiltonian, respectively, (c) INS spectrum as calcu- 
lated numerically from the microscopic Hamiltonian and the 
single-spin approach for a temperature T — 3 K, momentum 
transfer QR = 4, and Gaussian lines with a FWHM of 20 /xeV. 
The curves were scaled so as to coincide for the 2 — > 3 peak at 
0.3 meV. The inset shows the energy spectrum of the dimer. 



rameters, the energy spectrum consists of well separated 
spin multiplets as shown in the inset of Fig.^c). 

The second-order effective spin Hamiltonian of the 
5 = 3 ground-state multiplet is -ff = D[S'^ — 
1/3S{S + 1)] -I- B^Ol [as obtained from eq. ^ for 
a uniaxial cluster]. For the calculation of the pa- 
the first-order projection coef- 



rameters D and B^, 



ficient r2(l) and the second-order projection coeffi- 
cients r2,2(l) and r2.4(l), which enter D, Di and 
[eqs. (|23f) and (|24|) ]. need to be evaluated, or the two 
reduced matrix elements (5/2, 1/2, 3|lr2(5'i)||5/2, 1/2,3) 
and (5/2,1/2,3||T2(5^)||5/2,1/2,2), respectively [see 
eqs. (|15ll and H17() ]. Here, we have written the zero-order 
dimer states as |5'i5'25'M). The respective coefficients 
for ^2 are zero as a spin- 1/2 cannot exhibit a zero- field- 
splitting. Using standard ITO methods, the calculation 
is easy enough and one obtains D = 2/3I?i + A/GSDf/ A 
= 0.6709 K and = l/SlbDj/A = 0.2116 mK [the 
other parameters in eq. (|23|) are obviously zero] . If one 
fits the exact eigenvalues of the states of the S = 3 mul- 
tiplet to the effective spin Hamiltonian H, one obtains D 
= 0.6707 K and B^ = 0.1770 mK, showing that the ana- 
lytic second-order results give the corrections D — 2/3Di 
and i?2 correctly up to 5% and 16%, respectively. 

Next, the INS intensity, eq. (|3nj), will be evaluated. 
The parameters Di and D3, respectively, have been al- 
ready obtained in the above calculation of the effective 
spin Hamiltonian. For the further evaluation it is useful 
to note the special properties of a dimer, namely that 

ri(i) + ri(2) = 1, ri(i)/ri(2) - S1/S2, and ri,fe(2) - 

— ri^fe(l). The first two relations determine ri(l) — 5/6 
and 'ri(2) = 1/6, which yields f{Q) = 1 - 5/18J^(g) 
with J^{Q) = 1_- joiQR) + l/2j2{QR), and the third 
relation yields fk{Q) = 2/3Ti^kU)^{Q)- One is thus 
left with the calculation of the reduced matrix element 
(5/2,l/2,3||fi('5'i)||5/2,l/2,2) entering ri,fc(l), which 
is again easy enough. One finally obtains the result 
fniQ) = 1 - [5/18 - 4/27Dl/A{AP - l)]J'{Q), which 
is plotted in Fig. ^c) for the three possible transitions 
within the S = 3 multiplet. For comparison, the result 
of an exact numerical calculation is shown in Fig. ^b). 
The agreement of the second-order and exact results is 
good. In particular, the correct trend and order of mag- 
nitude is recovered, though the second-order results for 
fi{Q), which corresponds to the transition 1 — > 2, lie 
significantly too low at larger Q values [for QR ^ 1 the 
fi{Q) curve in Fig. ^a) is visibly shifted towards the 
fo{Q) curve as compared to Fig.QIb)]. 

These results demonstrate the considerations of the 
previous subsection, in particular the dependence on M 
is apparent. In the case of a dimer, the first- and second- 
order contributions exhibit identical Q dependencies, i.e., 
f{Q) and fk{Q) in eq. H30I) are both proportional to the 
same function !F{Q), since the spatial symmetry of a 
dimer completely fixes the functional dependence on 
In the general case, however, f{Q) and fk{Q) will exhibit 
different Q dependencies because of the different topolo- 
gies of the spin transition maps involved in f{Q) and 
fk{Q), respectively. 

The dependence of the interference factors on M ob- 
viously leads to changes in the intensity pattern of the 
transitions as compared to that calculated in the single- 
spin approach, or in first order, respectively (both the 
single-spin approach and the first-order result disregard 
the M dependence and are similar in this regard). This 
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is shown in Fig.^c), which presents the neutron energy- 
loss spectrum as calculated exactly and in the single-spin 
approach (the results were obtained numerically). The 
calculated intensities were normalized so as to coincide 
for the 2^3 transition at 0.3 meV. Clearly, the many- 
spin effects lead to a significant effect on the INS inten- 
sity, which is on the order of 10% for the present system. 
This is encouraging, as it suggests that such effects in- 
deed could be observable. It is also noteworthy that the 
many-spin effects on the INS intensity are much stronger 
than on the zero-field-splitting parameter I?, which is af- 
fected by only 0.6%. 

The above example allows to elucidate another impor- 
tant point. If the spin-5/2 center would additionally ex- 
hibit a single-ion forth-order anisotropy, which would add 
a term BiOl{Si) to the microscopic spin Hamiltonian, 
the fourth-order term in the effective spin Hamiltonian 
would become = l/315I?f/A + l/Wi. Clearly, an 
analysis of the energy spectrum or transition energies 
alone cannot distinguish between the second- and first- 
order contributions to coming from Di and Bi, re- 
spectively. In contrast, the many-spin effects on the INS 
intensity are governed by the spin mixing due to the Di 
term, while the Bi term is ineffective (see the discussion 
at the end of section [3 A). In other words, in the INS 
spectrum, Fig. ^c), the transition energies would shift 
with an additional Bi term, but the intensities would 
not be affected. Thus, a simultaneous analysis of energy 
spectrum and transition intensities allows, in principle, to 
separate the parts in the cluster fourth-order term which 
come from spin mixing and from single-ion anisotropy. 
This has great potential. With the advent of the single- 
molecule magnets the understanding of the origin of mag- 
netic anisotropy, and its eventual controlling, has grown 
enormously in importance. Fourth-order terms in the ef- 
fective spin Hamiltonian such as a Bf term are of partic- 
ular relevance because of their crucial role in the process 
of quantum tunneling of the magnetizationi^iSS INS stud- 
ies as the one just suggested have the potential to evolve 
into a powerful tool in this regard. 

As a second example, the so called Mn-[3 x 3] grid 
will be considered. This system is of some relevance for 
the present work as it provided a clear-cut demonstration 
of the importance of spin mixing via the observation of 
magneto-oscillations in the torque signalJfi In the Mn- 
[3 X 3] grid molecule, nine spin-5/2 Mn(II) ions occupy 
the positions of a regular 3x3 matrix, held in place by a 
lattice of organic ligands. The parameters of the micro- 
scopic spin Hamiltonian were determined from magneti- 
zation, high-field torque, and INS studies. ^O'^^"*-^^ The Mn 
ions within a molecule experience an antiferromagnetic 
interaction leading to a 5* = 5/2 ground state, which ex- 
hibits an easy-axis zero-field-splitting with D = —0.48 K. 
In the following, the two transitions 1/2 3/2 and 
3/2 5/2 within the S = 5/2 ground-state multiplet 
are the subject of interest. They have been observed by 
INS, though with poor statistics. 

Using sparse-matrix techniques, the interference fac- 



tors /i/2(Q) and /3/2(Q) characterizing the two transi- 
tions under consideration as well as the neutron energy- 
loss spectrum were calculated numerically from the full 
microscopic spin Hamiltonian [for the Hamiltonian see 
eq. (1) of Ref . Hol and the parameters given there]. The 
results are presented in Fig. |21 They once more demon- 
strate all the effects of spin mixing identified before. In 
particular, the interference factors shown in Fig.|2Ia) are 
again clearly different for the two transitions, i.e., are 
dependent on M. The exact intensities in fact differ by 
as much as 9% as compared to the expectation from the 
single-spin approach. The peak intensities in the neutron 
scattering spectra are accordingly affected too, as demon- 
strated in Fig. |2Ib). This example thus shows, that the 
effects of spin mixing on the INS intensities indeed can 
be of significant size in experimentally relevant systems. 
First hints for the relevance of such effects have been in 
fact found in a previous INS studyi^ 

The example of the Mn-[3 x 3] grid also allows to 
demonstrate the other effect of spin mixing, which was 
not displayed by the above dimer because of its partic- 
ular symmetry properties. In Mn-[3 x 3], the Q depen- 
dencies of the two transitions differ not only in magni- 
tude, but also exhibit different functional dependencies. 
This is evidenced in the inset of Fig. [2Ia), which shows 
the behavior of the Q dependence of the two transitions 
1/2 — > 3/2 and 3/2 5/2 at low momentum transfer 
in greater detail. A crossing of the two curves is evi- 
dent, which indicates that the functions f{Q) and fk{Q) 
in eq. (|30|) exhibit different functional dependencies on Q 
because of the different spin transition maps involved in 
/(Q) and h{Q). 

As a final example, the famous Mni2 cluster will be 
consideredi^ This single-molecule magnet consists of a 
tetrahedral core of four spin-3/2 Mn(IV) ions surrounded 
by a crown of eight spin-2 Mn(III) ions, and is char- 
acterized by a 5 = 10 ground state with a large easy- 
axis zero- field-splitting [D = -0.66 K, B^ = -34 /iK, 
Bj = ±43 /iK).22'29 Furthermore, starting at about 40 K 
a number of higher lying spin multiplets with 5 = 9, 10, 
and 11 exist, which are so low in energy that the lowest 
S ^ 9 multiplet actually overlaps with states of the zero- 
field-split 5 = 10 ground-state multiplet [see Fig. 7(a) in 
Ref. i57j . Because of this, one expects pronounced spin- 
mixing effects in this cluster. 

Using again sparse-matrix techniques, the energies of 
the states with Af = 0, . . . , 10 of the 5 = 10 ground-state 
multiplet, as well as the interference factors fo{Q) to 
hiQ) for the respective transitions between them, were 
calculated numerically from the microscopic spin Hamil- 
tonian. The Hamiltonian consisted of the isotropic ex- 
change terms appended by uniaxial single-ion anisotropy 
terms. The exchange part and the respective coupling 
constants were chosen as in Ref. ^3- Since the Mn(III) 
Jahn- Teller ions are expected to be the main source of 
magnetic anisotropy, the single-ion anisotropy was mod- 
elled by a term Dj^n{iii)^Sf ^, where the sum runs 
over the eight Mn(III) ions. From the first-order pro- 



13 



f„, (1/2 ^3/2) 

' - 8.9% _ f^^^ (3/2 ^ 5/2) 




0.10 0.15 0.20 

E (meV) 



FIG. 2: Results of numerical calculations for the Mn-[3 x 3] 
cluster, (a) Q dependence of the interference factors /iw for 
the two transitions 1/2 — > 3/2 and 3/2 —> 5/2 within the 
5=5/2 ground state multiplet as calculated from the micro- 
scopic spin Hamiltonian. The inset shows the low-Q range 
with greater detail, (b) INS spectrum as calculated from 
the microscopic Hamiltonian and the single-spin approach for 
a temperature T = 3 K, momentum transfer Q = 1.5 A, 
and Gaussian lines with a FWHM of 20 ^eV. The curves 
were scaled so as to coincide for the 1/2 —> 3/2 peak at 
0.097 meV. The inset shows a sketch of the zero-field-split 
5 = 5/2 ground-state multiplet of Mn-[3 x 3]. 



jection coefficients one estimates D — Q.201ADMn{iii)i 
and DMn(iii) = -3-3 K was chosen for the calculations. 
The obtained interference factors are plotted in Fig. 
The ten curves are basically indistinguishable and co- 
incide with the first-order contribution /(Q) [which is 
easy to calculate analytically from eq. (|29|) thanks to the 
tables given in Ref. l57l |. Thus, surprisingly, spin mix- 
ing is very weak in Mni2. This finding is further sup- 
ported by fitting the calculated energies of the states of 
the S — 10 ground state to the effective spin Hamilto- 
nian H = D[Sl - lliS{S + 1)] + Biol, which yields D 
— -0.653 K and Bl — -1.97 jjK {Bf is obviously zero in 
our model), i.e., yields a value for Bl which is one order 
of magnitude smaller than the experimental one (but at 
least of the correct sign). 

The absence of spin mixing in Mni2 is rather surpris- 
ing since from the overlapping spin multiplets one might 
have expected strong mixing effects. The explanation is 
found by realizing that transition matrix elements are 
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FIG. 3: Numerically calculated Q dependence of the inter- 
ference factors /o to /g for the ten transitions within the S 
= 10 ground state of the Mni2 cluster (for details see text). 
The curves essentially fall onto each other. The inset provides 
a closer view on the peak at about Q = 1 A. The individual 
curves can be distinguished, but the differences are very small. 

not only subject to the selection rules arising from the 
spin rotational symmetry (such as jS* — S"| = 0, 1, 2 for 
the spin mixing due to a bilinear anisotropy) , but also to 
those arising from the spatial, spin permutational prop- 
erties of the microscopic spin Hamiltonian. '^'^•'^•^ In fact, 
for the Mni2 cluster, or more precisely, for the micro- 
scopic spin Hamiltonian assumed for it,^^ it turns out 
that the two lowest lying excited spin multiplets have a 
different spatial symmetry than the S" = 10 ground-state 
multiplet, and thus cannot mix with it [so that the quan- 
tities Z)fe, which appear in both the effective Hamilto- 
nian, eq. it^ . and the effective spin operators, eq. l(?7|l . 
are zero for these states]. Obviously, a high molecular 
symmetry can effectively switch off a mixing with many 
of the higher lying spin multiplets. This teaches the im- 
portant lesson that the spatial symmetry (or coupling 
topology) of the cluster plays an important role in the 
spin mixing - and in all the effects related to it. With 
respect to the magnetic anisotropy of the cluster, i.e., the 
effective spin Hamiltonian, this implies that a deviation 
of the actual molecule from an assumed high symmetry 
can lead to substantially different values of the param- 
eters. This effect should be suspected as an important 
origin for higher-order terms in the effective spin Hamil- 
tonian (and for Mni2 orthorhombic distortions in fact 
were demonstrated recently, see Ref. ^2^. and references 
cited therein). Though significant spin mixing is absent 
in the assumed model of Mni2, it was discussed here as 
an example as it allows these important and general con- 
clusions concerning the role of spin mixing and the origin 
of magnetic anisotropy. 

VI. ELECTRON PARAMAGNETIC 
RESONANCE 

In the previous section, and in eq. H23|l of section Hvl 
we have considered a system in zero magnetic field under 



14 



rather general conditions. In this section, we additionally 
take into account the Zeeman term in the microscopic 
Hamiltonian, eqs. |2l or (j2J. Otherwise, we make the 
same assumptions as before, i.e., we assume Hi — J^ij 
Aij • Sj + ^ ■ Si • gi • B with traceless and diagonal 
Aij (we will refer to the two terms in Hi as Hd and Hz, 
respectively) . 

First, we discuss the effective spin Hamiltonian. Equa- 
tion H21|l of section IIVI can be applied directly. The 
only new point which needs to be addressed here is that 
Hi now consists of two terms instead of one. Accord- 
ingly, the first-order part of the effective Hamiltonian 
H also consists of two terms, namely S • D • S and 
/isS • g • B (these two terms will be denoted symboli- 
cally as SDS and SgB), where D = J^ij ^2{'ij)-Aij and 
g — X]iri(i)gi. For the second-order part, however, 

the four terms PHdP'HdP, PHzP'HzP, PHdP'HzP, 
and PHzP' HjjP have to be evaluated. 

The term PHdP' H^P has been calculated already in 
section IIVI eq. (|23|l . It renormalizes the bilinear zero- 
field-splitting term (the symbol SDS shall be understood 
as to include this effect) and produces fourth-order terms 
(abbreviated symbolically as S^). 

The term PHzP' HzP, which is quadratic in the mag- 
netic field, adds a fc = and k = 2 contribution [and a 
k = 1 contribution, which, however, is related to vector 
products bj x hj and thus smaller by a factor (Sg/go)^; 
it is neglected here]. The k — contribution produces 
what is known as Van-Vleck or temperature-independent 
paramagnetism (TIP). It may be written as —^xoP^ (or 
symbolically as B^). The TIP parameter xo can be cal- 
culated straightforwardly from eq. H21|) . The k = 2 con- 
tribution can be put in the form S • A • S, with the ten- 
sor A composed of the (reduced) magnetic fields only: 
A = ^j^- hi ■ hJ . Thus, this contribution adds to the 
bilinear anisotropy, but the anisotropy depends on the 
field direction and magnitude (it is denoted symbolically 
as SB^S) 

The terms PHdP'HzP and PHzP'HdP, which both 
are linear in the magnetic field, produce a A: = 1 and 
fc = 3 contribution (the k = 2 term cancels out as wave 
functions may be chosen as real). The k = 1 contribu- 
tion, since it is linear in both the magnetic field and the 
spin operators, adds a term ^bS • Ag • B to the effective 
Hamiltonian [where Ag can be obtained again straight- 
forwardly from eq. (|21|l ]. Thus, this part modifies the 
effective g tensor, making it more anisotropic (the sym- 
bol SgB shall be understood as to include this effect). 
The fc = 3 contribution, however, is unusual as it is of 
third order in the spin operators. It cannot be put into a 
bilinear form as the other terms (it is symbolically writ- 
ten as S^B). 

Altogether, the second-order effective spin Hamilto- 
nian consists of the six terms SDS, S^, SB'^S, S^B, 
SgB, and B^. The above considerations seem pretty 
obvious, but appear to have been disregarded in all ac- 
tual analyses of EPR results on molecular spin clusters. 



The effective spin Hamiltonian used there consisted of 
the zero-field-splitting terms (up to any desired order) 
and the Zeeman term (with a possibly anisotropic g ten- 
sor). That is, in the present context only the terms SDS, 
S'^, and SgB have been considered. The B^ term does 
not affect energy differences and thus can be omitted in 
spectroscopic work (though it has been demonstrated to 
be of importance for a correct interpretation of high-field 
thermodynamic datai^). However, we do not see a justi- 
fication for omitting the remaining field-dependent terms 
SB^S and S^B, especially in high-field EPR (HFEPR). 

In low-field EPR {X- and Q-band EPR), the terms 
SB^S and S^B might affect the analysis little, at least 
one can hope so. But in HFEPR, which is actually mo- 
tivated by the aim to work in the limit Hz ^ Ho, the 
effects introduced by these two terms should be signif- 
icant. This might be objected to with the general rea- 
soning that these terms, as they are of second order in 
the perturbation theory, are reduced by 1/A and thus 
always small. However, whenever there are reasons to 
expect that e.g. the parameters B^ are significantly in- 
fluenced by spin mixing (as it is often the case and has 
been demonstrated explicitly for e.g. the single-molecule 
magnets Feg, Fe4, and the molecular wheel CrgiiiSiS), 
one also should expect significant effects from the terms 
SB^S and S^B. 

A more detailed analysis of the resulting EPR spectra 
would be beyond the scope of this work even for the sim- 
plest situations. We want to close with some comments. 

From the above discussion it should be clear that with 
the current practice of neglecting non-linear field depen- 
dencies in the spin-Hamiltonian, the parameter values 
determined by HFEPR and INS may differ substantially 
(for the system Fes, e.g., it has been reported that the ax- 
ial parameter D as determined by INS is 6% larger than 
that obtained by HFEPR^^). However, the non-Zeeman 
like field dependence of the operators, in turn could be 
developed into a powerful tool to analyze the many-spin 
effects: a careful comparison of spectra recorded at differ- 
ent fields/frequencies provides direct information about 
the second-order field-dependent terms in the effective 
spin Hamiltonian, and thereby about the effects due to 
spin mixing. 

The calculation of the effective spin operators Sia is 
simpler than the calculation of H in the sense, that the 

second-order contributions to Sia can be calculated in- 
dependently for each term in Hi, and then be added up. 
Thus, the results of the previous section, eq. H27|) . pro- 
vide directly the second-order contributions due to Ho- 
The calculation of the second-order contribution of the 
Zeeman term is, by following the same line as in the pre- 
vious section, easy enough, but the result is again pretty 
lengthy and therefore not reported here. 

As for INS, a general result is possible for the EPR in- 
tensity if one considers only the first-order contributions 
in the effective spin operators. In linear- response theory. 
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the absorbed EPR power is given by 
P(B) = du:, T) ^ ^P_(B)J(c. - £2!i_±!i), 

nm 

(31) 

and 

P„™(B) - ^^C6;.^(n|^,„|m)(m|^,^|n), (32) 

ij a/3 

with C(a;,r) = - e-''^")/(87rS), and the (reduced) 
microwave field = /isgi • B' (which is to be contrasted 
to the static magnetic field B, or b^). Insertion of Sia = 
Ti{i)Sci leads directly to the first-order result 

P„i^(B) ^Y.W^\SA^)^MSp\n), (33) 

a/3 

with the field b' — /i^g • B' [we note that for the case 
gi =51, the EPR intensity becomes proportional to the 
INS intensity at Q = 0, see eq. (|^ ]. Thus, as expected 
and mentioned at the end of section m the many-spin 
nature of the wave functions does not affect the relative 
intensities of the EPR transitions in first-order perturba- 
tion theory. Up to this order, the single-spin approach 
is justified. But as in the case of INS, second-order ef- 
fects due to spin mixing can be significant too. How- 
ever, in view of the inaccuracies introduced already by 
the use of an incomplete spin Hamiltonian in the analy- 
sis of (IIF)EPR results, see the above discussion, we are 
not pursuing this issue further here. 

VII. CONCLUSIONS 

In this work, we have investigated the influence of spin 
mixing, i.e., a mixing between the different spin multi- 
plets in the spectrum of a Heisenberg spin system, on the 
spectroscopic intensities of transitions within a particu- 
lar spin multiplet. This situation is encountered in many 
molecular nanomagnets, with systems such as the single- 
molecule magnets as canonical examples. The spectro- 
scopic intensities were calculated perturbationally, with 
the isotropic Heisenberg exchange terms as the unper- 
turbed part. Spin mixing arises then from the anisotropic 
terms in the microscopic spin Hamiltonian, and its ef- 
fects were calculated up to second-order in perturbation 
theory. The consequent use of the spherical tensor and 
ITO formalism yielded a number of useful, and easily 
evaluated, equations. The resulting mixing effects were 
discussed for the specific examples of INS and EPR, with 
strong emphasis on the former technique. 

As a main result it was observed that the transition in- 
tensities are affected individually for each transition, as 
compared to the results of a first-order calculation or the 
single-spin approach, respectively, which has been the 
approach used in all spectroscopic studies on the split- 
ting in spin ground states so far. The transition intensi- 
ties in fact were found to provide additional information. 



which in some cases may not be obtained otherwise. For 
instance - and this is probably the most important ex- 
ample for that - we showed that the careful inclusion 
of the transition intensities in the analysis of the ob- 
served spectra allows to differentiate between the two 
origins of higher-order magnetic anisotropy terms in the 
ground states, namely spin mixing or single-ion higher- 
order terms. Such studies would thus allow an unprece- 
dented insight into the microscopic origin of magnetic 
anisotropy in spin clusters, which ultimately could help 
to achieve a tailored design of magnetic molecules. We 
also showed, that the symmetry of the cluster, or its 
absence, can have a profound impact on the resulting 
anisotropy of the spin clusters. 

As a summary, in this work we established a complete 
toolbox for the analytical analysis of the spin-mixing ef- 
fects, including a general equation for the construction 
of the second-order effective spin Hamiltonian, and ex- 
plored the various resulting effects, as well as the kind of 
information which may be obtained from them. 

Finally, we should address the experimental feasibility 
of such investigations. The considered examples showed 
that changes in the intensities due to spin mixing may be 
on the order of 10% or so. With the INS spectrometers 
available nowadays, such effects should be detectable for 
magnetically strongly scattering species, which would be 
typically low-nuclearity spin clusters like dimers. Such 
systems in fact would be excellent model compounds for 
a demonstration of the proposed mixing effects or de- 
velopment of the methodology. For the most interest- 
ing molecules, however, like the single-molecule magnets, 
such studies seem to be prohibitive, if not impossible, at 
the moment. For EPR the situation is rather similar. 
However, considering the enormous improvements of the 
INS and EPR techniques in the last decade, the outlined 
studies should become possible in the future. We hope, 
that this work will stimulate development efforts in this 
direction - the new opportunities which would be gained, 
in our opinion, would more than reward these efforts. 



APPENDIX A 

Some comments shall be added concerning the Carte- 
sian and spherical representations of a spin Hamilto- 
nian. The general form of each term is predetermined 
by the fact that a Hamiltonian has to transform as a 
scalar under rotations in space. The scalar product of 
two spherical tensors has the gener alform E,t^.''^*#', 
or Eg(-l)''C^-9VJ''^ since ult'^* = {^lyu'^^l (here U 
and V are arbitrary spherical tensors, e.g., ITOs). The 
reader may notice this motive in many equations in this 
manuscript. In the ITO representation,'*^ the transfor- 
mation properties thus enforce the spin Hamiltonian to 
appear - in the notation of this work - as a sum of terms 
\Ttke^^{—lYT'^^Tq^\s), i.e., asasum of spherical scalar 
products. 
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For instance, for a bilinear term follows 

s,.A.s,= J2 E^f'^*w''' (Ai) 

fe=0,r,2 q 

The association of the fc = 0, 1, 2 contributions in this ex- 
pression with the isotropic, anti-symmetric, and symmet- 
ric&traceless parts of the (Cartesian) tensor A is stan- 
dard. The generalization of this concept to the other 
terms in the spin Hamiltonian is obvious. For instance, 
the Zeeman term fis^i ■ gi • B, which can be written as 
Si ■ hi with hi = ^Bgi • B, has to become the scalar 
product of the spherical tensor Tt\hi) and the tensor 
operator ft\i), both of rank one. 

Any arbitrary operator 0{S), which is a function of 
spin operators, can be expanded into ITOs according to 
^(5') = Y^i^qCkqfq''\s). Section lTlI Bl provides a general 
procedure to find the coefficients Ckq- For the bilinear 
term, eq. (|A1|) . one would find Ckq = (— I)tI^^(A) in 
full agreement with the result concluded from the trans- 
formation properties. 

The bilinear term, eq. IjAip . shall be discussed fur- 
ther. Let us consider a diagonal and traceless tensor A 
[that is, the fc = 0, 1 contributions in eq. (|A1() are zero]. 
The components of the spherical tensor Tq ' (A) associ- 
ated with A are T^\a) ~ {2Kzz — ^xx — ^yy)/V^ and 
T^^\k) = T^^iiA) = {Axx - Ayy)/2. The g = ±1 com- 
poncnts vanish. This suggests to define the parameters 

D = {2A,,~AxX~Ayy)/2, 

E = (A,, -A^^)/2, (A2) 

so that T^^\A) = ^j2/lD and T^a'C^) = ^ order to 
parameterize the interaction described by Si • A • S j . 
In Refs. 1^3 and 0, A was parameterized by J" — 

(2A^^ - Axx - Ayy)/V6 and J" = [Axx ~ ^yy)l^- 

J" agrees with T(j^^(A), but J" differs from Tfl{A) 
by a factor of \f2. We prefer our choice of parame- 
ters because of two reasons. Firstly, it is consistent with 
the usually chosen parametrization A^^. = — 1/3-D + £', 
t^yy = -1/3D - E, and A^^ = ^jZD {which for i = j 

leads to the spin operators D\Sl^ — \Si{Si -1-1)] and 
E{Sf^ — Sf )}. Secondly, it ensures the above compact 



form eq. (jAl|) of the Hamiltonian, which simplifies calcu- 
lations enormously as the rotational invariance properties 
are inherently manifest. 



APPENDIX B 

Table m lists the most often needed ITOs explicitly 
expressed in operator form. With these expressions and 
the Wigner 3j symbol 

(S S fc\ / (25)!/(25-fc)7~ 

1^5 -5 0^ V (2^-f l + fc)!/(2^)! ^ ' 



TABLE II: Some irreducible tensor operators which are often 
needed in practice. O4, O4, and Of are Stevens operators^- 
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it is straight forward to calculate the reduced matrix 
elements (S-Hf (3)(5)||5) and (5||f W(5')||S') not usu- 
ally found in textbooks. For fc < 4 one obtains 

{S\\f(''HS)\\S) ^ V^Hiif with to = 1, ii = 4, 
t2 = 24, ^3 = 160, and ^4 = 1120. 



APPENDIX C 

For calculations an explicit expression of the INS for- 
mula eq. H26|l is often needed, which is given here for con- 
venience. With /„™(Q) = E^, ^r(Q)^j (Q)^nm,.j(Q), 
eq. reads 



Inm.ijiQ) — T:jo{QRij){SixSjx -\' SiySjy -\- SizSjz) ~\' —j2{QRij) {'^SizSjz SixSjx ^iy^jy) 

o D Hi 



1 R —R R R ~ ~ 

^h{QRij) —{SixSjx — SiySjy) + j2{QRij) {SixSjy + SiySjx) 



. Rij ,xRij ,z / Pi , Pf pi \ , ■ / ^ 7-i \Rij.vRi 



+j2{QRij) '' p2 {SixSjz + SizSjx) + j2{QRij) — p2 ^"''^ (SiySjz + SizSjy). (CI) 



^ij ^ij 



r 



For real wave functions, which is the situation mostly r.h.s. vanish and can be dropped, 
encountered in praxis, the fourth and sixth term on the 
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It is actually often useful to have this formula in an- 
other form. Since the -dependent prefactors do not 
change upon the substitution ij — > ji, one may replace 
each SiaSjfj by {SiaSjp + SjaSip)/2. If one further ex- 
presses the spin operators Six and Siy by the ladder op- 



erators Sf, uses {n\Sf\m) = {m\S^'\n)*, and introduces 
the abbreviations i+ = {n\S^\m) and similar for the z 
and — components, one obtains the handy equation 



Inm,tjiQ) = -joiQRij)Re[2i^jl + i+jl + + —j2iQRtj) — '-^Re[4i,jl - i+j^ - 

+J2(Q%) ^^^"^^ Re[(z+ + ^-)J:] + 32m^,) \f'''' MiH (C2) 



where Re[] and Im[] denote the real and imaginary 
parts, respectively. This form of eq. (PHI is also use- 
ful if the spin matrix elements are calculated by the 
ITO techniquaSii since the corresponding ITO expressions 
iz = (n\f^'^\i)\m) and i± = =F\/2(n|f^^\i)|m) are read- 
ily inserted. 
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